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ABSTRACT 

f^ ■ In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far- 

r\l IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project 

, GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at 

I— 1^ several key gas-lines, including [01] 63.2 and 145.5yum, [Cll] 157.7yum, CO 72.8 and 90.2yum, and 0-H2O 78.7 and 179.5yum. We only detect the 

Q^ [01] 63.2yumline in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and '^''^CO data 

ry^ obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We 

. model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo- 

yS^ chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous 

CIh structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as 

the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results 

are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 ± 2.0 x 10"^ Mq is still present in this disc, in agreement with 

earlier CO observations. 

Cu Key words. Stars: pre-main sequence; (Stars:) planetary systems; (Stars:) circumstellar matter; protoplanetary discs; Infrared: planetary systems 
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>■ 1 . Introduction to a circumstellar (CS) disc. Submillimeter Array (SMA) obser- 

00 vations show a disc in Keplerian rotation with radius r = 235 AU 

(M Giant gas-planets form in protoplanetary discs within the first ^^^ inclination i ~ 13°(Raman et al. l2006l l and a total gas mass 

00 10 ]yiyr after protostar formation. Therefore the amount of gas ^f 0.6-3.O xlO^ M^ (Panic et al. "2008^. For a Herbig Ae, the 

rn ; that is present in a disc at a given time is very important, as it ^t^. j^ ^^^^^^1 ^^ j^ ^as a small near-IR excess (e.g. Dominik et 

V^ determines whether these planets can still be formed. The study ^j 2003 ). Furthermore, SWS/ISO spectra revealed that the sili- 

O of the dust IS equally important, as it witnesses the first steps of ^ate 1 fim feature, detected in the majority of the He rbig Ae/Be 

O planet formation. A lot is known about the dust in protoplane- ^^^^ (^^jy g q^^ ^f 53 j^^j^ j^e feature, Juhas z et a Usubmittedi 

^H . tary discs thanks to the Space Observatories ISO and Spitzer: j^ j^^j^j^g ^^ h^, 159142 (Meeus et aL lIOOTT l. This absence can 

L; ; different degrees of dust processing were observed with no clear ^g explained if the silicate grains are either too large or too cold 

ITj '"^'''^'°" ^° '^^"'''" properties, while grain growth could be related ^^ gj^it at 10 ^im (Meeus et al.l2002l. On the other hand, features 

P^ to the disc structure (e.g. Sicdia-Aguilar et al. 2007, Meeus et ^f polycyclic aromatic hydrocarbons (PAHs), which can be ex- 

H ; al. 12009 ). Gas is more difficult to observe, as the special lines ^jted ^y uv photons, were clearly detected with ISO and Spitzer. 

. .^ , are not very strong. Most gas studies are based on CO lines in 3^^^^ on the near to far-IR ratio, Grady et al. (I2007D suggested 

the near-IR (warm gas) or in the mm, where CO freeze-out onto ^^^^ j^e inner region has already cleared some material, and that 

grain s is a complicating factor (e.g. Bnttain et al. ^2007] Dent t^g ^^^^^j. ^nd outer disc are not coupled. However, this might 

et al. 120051. However, as gas lines are stronger in the far-IR, ^1^^ ^e explained by a reduced opacity due to grain growth. 

we expect this field to experience major breakthroughs in the j„ ^^^-^ ^^ ^1^^^ ^j^^ ^^^^ Herschel observations of 

coming years based on data obtained with the Herschel Space j^j^ ^^^g^^^, and use radiative nansfer and chemical models to 

Observatory (Pilbratt et al. 2010), which provides sensitive far- ^^^^^^^^ j^g g^, properties of the outer disc. 

IR photometry and spectroscopy. 

The source HD 169142 is a young - age 6^3 Myr (Grady et al. 

I2007I I - Herbig Ae star with an IR to millimetre excess attributed 2 Herschel PACS Observations 



* Herschel is an ESA space observatory with science instruments Our observations are part of the Herschel Open Time Key 
provided by European-led Principal Investigator consortia and with im- Project GASPS (P.I. Dent, see Mathews and Thi, this iss ue, and 
portant participation from NASA. Dent et al. 120101 ). We obtained PACS (Poglitsch et al. 120101 ) 
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Fig. 1. Observed [OI] 63yum emission line (black) together with 
our line fit (red; parameters listed in TablelTJ. 

Table 1. Summary of the PACS targeted lines and the CO lines 
that were observed with the SMA (Raman et al. 120061 Panic et 
al. 120081) . We list the central wavelength, continuum and line 
flux or upper limits (three sigma) in the case of a non-detection. 
Between brackets we give the one sigma statistical error for the 
detections. 



Species 


Transition 


A, 


Continuum 


Line Flux 






(fim) 


(Jy) 


(10-'** W/m^) 


[OI] 


'Pi ^ 'P2 


63.18 


18.89(0.13) 


71.7(3.8) 


[OI] 


'Po ^ 'Pi 


145.53 


13.76 (0.03) 


< 10.4 


[CII] 


P3/2 -^ "Pl/2 


157.74 


14.48 (0.03) 


<6.4 


0-HoO 


2l2 -» loi 


179.53 


13.12(0.05) 


<8.8 


o-H.O 


423 -> 3,2 


78.74 


17.82 (0.32) 


< 10.6 


CO 


36^35 


72.84 


17.92 (0.04) 


< 15.5 


CO 


29^28 


90.16 


17.00 (0.05) 


< 10.6 


l^CO 


2^ 1 


1300.40 


0.169(0.005) 


0.093 (0.004) 


"CO 


2^ 1 


1360.22 


0.169(0.005) 


0.048 (0.004) 



photometry (obsid 1342183656, Point Source, 400 s) and spec- 
troscopy (obsid 1342186309, PacsRangeSpec, 5150s and ob- 
sid 1342186310, PacsLineSpec, 1669 s). The spectroscopic data 
were reduced with the developer build version 3.0.1212 of the 
Herschel Interactive Processing Environment (HIPE; Ott l2010| l. 
using standard tasks provided in HIPE. In order to conserve the 
best signal and not to introduce additional noise, we only ex- 
tracted the central spaxel and corrected for the flux loss with an 
aperture correction. Furthermore, we applied another correction 
factor as supplied by the PACS team (1.3 in the blue, and 1.1 in 
the red), to obtain a flux calibration accuracy of 40%. IVlore de- 
tails on the data reduction can be found in Mathews et al. (2010). 
In Fig. [T] we show our only detection, the fine structure line [OI] 
63.2 jum. We also searched for other gas lines that were targeted 
in our spectroscopy, but did not detect: [OI] 145.5 //m, [CII] 
157.7yum, CO 72.8 and 90.2yL(m, nor 0-H2O 78.7 and 179.5;um. 
In Table [1] we list the measured line properties. We note that 
in the LWS/ISO spectrum of HD 169142, the [CII] 157.7 yum 
line was detected, but LWS uses a larger beam (~ 40") than the 
size of the central spaxel of PACS (9'.'4 by 9'.'4). Therefore we 
also analysed the other 24 spaxels of the IFU and found a ten- 
tative detection in spaxels 13 and 14, which are to the West of 
HD 169142. The cause of this emission (if confirmed) is cur- 
rently unknown. 



Photometric observations were obtained in Point Source 
IVIode for both the blue (70 /im) and red (160/im) bands, and 
reduced using HIPE v2.3.4 with the Point Source Observations 
pipeline. The source is spatially unresolved, with FWHM 5" at 
70 /im and 11" at 160/^m. We derived photometry using an aper- 
ture of 21", and applied an aperture correction from the PACS 
PhotChopNod Release Note (Feb. 22, 2010). This gives a flux 
of 27.35 + 0.03 Jy at 70//m and 17.39 + 0.05 Jy at 160 //m, 
with a flux calibration uncertainty of 5% in the blue, and 10% 
in the red. Applying the flux uncertainty of 40% for the spec- 
troscopy, our photometry is consistent with the continuum fluxes 
measured from the spectral scans. Given the smaller error in the 
photometry, however, we give preference to those fluxes, and ex- 
pect to reduce the spectroscopic flux uncertainties in the future. 



3. Determination of stellar parameters 

Because previous studies list conflicting effective temperatures, 
we re-estimated the fundamental parameters. Photometry ex- 
cluding the contribution from the cool companion (at a distance 
of 9"3; Grady et al.l2007l ) was collected from Cutri et al. l2003l l. 
Zacharias et al. ( I2004I I and Sylvester et al. ( 1996). We also anal- 
ysed optical spectra, taken with CES/CAT with a resolution of 
~ 60000. Temperatures were selected by simultaneously fitting 
the lines from neutral and singly ionised species. The synthetic 
stellar spectra were computed using the ATLAS9 and SYNTHE 
codes by Kurucz (1993). x' tests show that models between 
Teff = 7500 K, \ogg, = 4.0, [Fe/H]=-0.50 and T^« = 7800 K, 
Iogg»=4.1, [Fe/H]=-0.25 match our data. 

Furthermore, the International Ultraviolet Explorer (lUE) 
obtained five spectra of HD 169142: one between 1200- 1900 A 
(SW) and four between 1900-3200A (LW). The SW spectrum is 
unusable for measurements below 1 650 A, but suggests emission 
lines of 01 (1304 A), CII (1335 A) and CIV (1550 A) - unfor- 
tunately, most are dominated by bad pixels. The LW spectra do 
not show any variability (over a period of 4.5 months), match the 
stellar photosphere, and do not show emission features. Overall, 
the lUE data (obtained from the INES INTA archive) do not give 
evidence for a UV excess. 



4. Analysis 

4. 1 . Dust and continuum modelling 

Far-IR lines emerging from a CS disc are affected by the stellar 
UV irradiation, disc mass and geometry, dust size and compo- 
sition, as well as PAH abundance. In order to interpret the lines 
observed with Herschel, it is crucial to first obtain a solid knowl- 
edge of the dust structure in the disc, based on as many observa- 
tions as possible, like SED, scattered light images and visibili- 
ties. Each of these observations provide complementary views of 
the disc structure and the dust properties. The disc model is cal- 
culated with the Monte Carlo radiative transfer code MCFOST 
(Pinte et al . 120061120091 1, as outlined below. 

We consider an axisymmetric, slightly flared density struc- 
ture with a Gaussian vertical profile, assuming power-laws for 
the surface density and scale height. We assume homogeneous 
and spherical dust grains (Mie theory), with sizes distributed ac- 
cording to the power-law f{a) oc a^^-^ between amin and dmax, 
where a is the grain radius. The dust is assumed to be well-mixed 
with the gas, i.e. the dust/gas ratio is constant throughout the 
disc. The star is reproduced by a uniformly radiating sphere with 
previously determined parameters Teff - 7800 K, log^» =4.1, 
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Table 2. Parameters of the best fitting dust model. 
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Fig. 2. Best fit with MCFOST to the observed data. Shown is 
the photometry obtained from the literature (black), Spitzer IRS 
spectrum (red), PACS photometric observations (blue crosses), 
and PACS continua derived from the spectroscopic observa- 
tions (red crosses). The green line represents the best MCFOST 
model, in red the atmosphere model. The blue dot-dashed and 
pink dashed lines are the contributions (scattered light and ther- 
mal emission) from the inner and outer disc, respectively. 



[Fe/H] - -0.25, and R, = 1.6 Rq. Parameters are adjusted to si- 
multaneously fit the SED, the Spitzer/IRS spectrum (Sloan et al. 
[2005), the 1 . 1 fim HST image (Grady et al. 2007) and the 1 .3 mm 
SMA visibilities (Panic et al. 120081 1. The scattered light image 
mainly constrains the flaring index to a low value (around 1.0) 
and the mm visibilities indicate a surface density varying as r"' . 

We reproduce the 1.3 mm emission with a dust mass of 
1.5 X 10""* Mq with a grain size distributiorQ between a^nm - 
0.03 jum and a^ax = 1 cm. The weak emission around lOfim im- 
plies that there is a discontinuity in the disc t-1 surface, possi- 
bly due to a shadowed area or a gap. In this paper, we explore a 
solution with a gap. For the dust composition, we use a mixture 
of 70 % silicates (Draine et al. .2003} and 30 % amorphous car- 
bon (ACAR sample; Zubko et al. 1 19961 1. and calculate the effec- 
tive optical index with the Bruggeman mixing rule. The strong 
PAH features around 10 /vm can be reproduced with a low PAH 
abundance due to the low continuum emission in that region. We 
modelled the emission with a single grain-size (54 carbon atoms, 
positively ionised) and assume that the PAHs are uniformly dis- 
tributed in the outer disc. We did not try to reproduce in detail 
the various bands observed with IRS, but rather constrain their 
abundance. We obtain MpAH/Mdust = 5 x 10""*, corresponding to 
/pah - 0.03 (where /pah is the uniform PAH abundance rela- 
tive to a standard ISM abundance of 10"^^^ PAH particles/H- 
nucleus, otpah = 667 amu, and gas/dust - 100). 

In Fig.|2]we show the best fitting model on top of the SED. 
Our model provides good constraints on the outer disc (mid-, far- 
IR and mm emission): dust mass, scale height, flaring index and 
surface density profile, as well as dust properties and amount 
of PAHs. However, the inner disc (near-IR emission) remains 
poorly constrained. We tried to improve the fit by moving the 
inner radius closer to the star {ri„ - 0.06 AU), but then the tem- 
perature is too high (2400 K) for dust, even carbon, to survive. 
Because we lack simultaneous visible and near-IR photometry 
as well as spatially resolved near-IR data, we did not further im- 



Parameter 


Inner Disc 


Outer Disc 


Hn (AU) 


0.1 


20 


'-OU, (AU) 


5 


235 


surface dens, exp e 


-1.0 


-1.0 


flaring exponent /? 


1.05 


1.00 


ref scale height ho 


0.07 AU @ 1 AU 


12.5 AU@ 100 AU 


Md„st (Mq) 


2 X 10-' 


1.5 X lO-t 


/pah 





0.03 



Table 3. Predicted line fluxes for different models, compared 
with the observed fluxes. The best fit is obtained by model#3. 
All models in this table have turbulent line broadening Vturb = 
0.15km/s. 





model #1 


model #2 


model #3 


observed 


gas/dust 


1 


100 


33 


- 


/pah 


0.01 


0.0055 


0.0087 


- 


/uv 


0.005 


0.0 


0.0 


- 


Line 




Line Fluxes [10-'*'W/m2] 




[01] 63.2fim 


154 


71.6 


71.6 


71.7 


[01] 145.5 ;um 


5.17 


10.1 


7.01 


<10.4 


[CII] 157.7/im 


4.58 


0.04 


0.06 


<6.4 


0-H2O 179.5 yum 


5.66 


5.15 


1.76 


<8.8 


'^CO 2 -^ 1 


0.060 


0.092 


0.093 


0.093 


"CO 2^ 1 


0.011 


0.059 


0.048 


0.048 


'^CO/'^CO 


5.69 


1.55 


1.92 


1.94 



' Note that for these size distribution parameters, only 31.5 % of the 
dust mass is in grains < 1 mm, the rest (68.5%) is in larger grains. 



prove the fit to the inner disc. However, as most of the mid-IR 
and mm lines originate from the outer disc, this is not critical for 
the present paper 

4.2. Gas and line modelling 

The best model reproducing the continuum observations is now 
fed into the gas thermo-chemical code ProDiMo (Woitke et 
al. I2OO9I 1 to calculate the chemical and gas temperature struc- 
ture in the disc and to predict the line fluxes, following the 
pipeline described in Woitke et al. (120101 see their Fig. 1). This 
final modelling step has additional free parameters such as the 
dust/gas ratio. We computed models for the following param- 
eters: gas/dust mass ratio e [1, 100], turbulent broadening with 
Vtuib e [0, 0.15] km/s, PAH abundance /pah e [0, 0.06] and stellar 
UV excess /uv - Luy/L-t, e [0, 0.005] (for details, see Woitke et 
al. 120 1 Oi l. Due to the stellar UV irradiation, the models gener- 
ally result in disc surface-layers where the gas is warmer than 
the dust. Table [3] shows the calculated line fluxes for a few se- 
lected models. An important result from our modelling effort is 
that the gas heating by PAHs plays a central role for the line 
flux predictions, in our case useful for the temperature sensitive 
[OI] 63.2 /im line. For each selection of the parameters dust/gas 
and /uv, we tuned the PAH abundance until a fit with the ob- 
served [OI] 63.2 fiin line flux was obtained, if possible. However, 
models with UV excess /uv =0.005 result in a much too hot gas 
and hence too strong gas emission lines, even for /pah - 0, and 
even if we decrease gas/dust — > 1, as in model #1. For the stan- 
dard ratio gas/dust - 100, we found model#2, which fits the 
[OI] 63.2 /im and '^CO line fluxes, but the '^^CO line is slightly 
too high (2.5 sigma). A better fit is obtained with gas/dust = 33 
(model #3), which simultaneously fits all three detected lines, 
and agrees with all other line upper limits. However, we empha- 
sise that apart from the uncertainties in the observations, system- 
atic uncertainties in the physical description of various chem- 
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Fig. 3. Spatial origin of the various gas emission-lines. The two 
vertical white dashed lines indicate 15% and 85% of the radially 
cumulative face-on line flux, respectively, i.e. 70% of the line 
flux originate from within the two white dashed lines. The red 
lines mark the cells that contribute most to the line flux in their 
vertical column. 



ical and radiative processes in the models render a proper gas 
mass determination difficult. The observational uncertainty in 
the '^CO line flux - which appears to be a dominant gas mass 
tracer - translates to a ratio range gas /dust - 22 to 50. An in- 
crease of the [OI] 63.2 //m flux by 40% (to account for the cali- 
bration uncertainty and match the photometry) would indicate a 
higher gas temperature. In our modelling, this is best interpreted 
by a larger PAH abundance (/pah = 0.02), as the other more rel- 
evant parameter, /uv, is constrained by the non-detection of the 
[CII] 157.7pm line. 

The spatial distribution of the '^^'^CO molecules predicted 
by our chemical model is consistent with the interferometric ob- 
servations by Panic et al. (I2008I I. Figure[3]shows a detailed anal- 
ysis of the spatial origin of the various emission lines, as de- 
rived from model #3. The [OI] 63.2 //m and 145.5 yum lines orig- 
inate from a radial disc region extending from the beginning of 
the outer disc at 20 AU to about 75 AU and 60 AU, respectively. 
Both lines are optically thick (rune ~ 10 - 100) and come from 
relative heights z/rx; 0.2 - 0.3. The [CII] 157.7 yum line is op- 
tically thin and extraordinarily weak in this model, due to the 
lack of stellar UV photons < llOnm that are capable of ion- 
ising carbon (/uv = 0). The '^CO and '^CO / = 2 ^ 1 lines 
probe the conditions in the outermost disc regions 43 - 185 AU 
and 52 - 190 AU, respectively, and are extremely optically thick 
(Tiine ~ 100 - 5000). Whereas the '^CO line comes from the 



PDR-like CO surface z/r x 0.4, the more transparent '^CO line 
probes slightly deeper layers z/r ^ 0.2. These findings are quite 
stable and robust among all calculated models. 

5. Conclusion 

We presented the first PACS observations of HD 169142 and 
showed the unique capability of Herschel to obtain an indepen- 
dent gas mass determination by simultaneously modelling the 
atomic fine-structure lines that are temperature sensitive (due 
to their high excitation energies), when combined with ground- 
based observations of '-'''^CO-lines. We showed that the obser- 
vations are consistent with a simple model of a disc hosting a 
gap, with a continuous structure between 20 and 200 AU in 
which the gas and dust are well-mixed, but stress that this so- 
lution is not unique. We determined the location of the emitting 
species, and constrained the gas/dust ratio to be ~ 22-50. We also 
derived that the UV excess emission, if present, must be lower 
than /uv = 0.005. We derive a gas mass between 3.0 and 7.0 x 
lO""* Mq, in agreement with the gas mass determination based 
on the CO lines alone (Panic et al. 2008). Despite indications 
that the HD 169142 disc is transitional, it still is a gas-rich disc. 
Furthermore, although there is dusty material close to the central 
star, the UV excess induced by accretion is extremely weak, if at 
all present. We will apply a similar modelling strategy for future 
GASPS observations, with the aim to study gas dissipation in a 
wide range of protoplanetary discs. 



Acknowledgements. We thank the PACS instrument team for their dedicated 
support, O. Panic for a discussion on her data and A. Juhasz for providing the 
SpitzerflRS data. C. Eiroa, G. Meeus, I. Mendigutia and B. Montesinos are pailly 
supported by AYA 2008-01727; J.M. Alacid and E. Solano by AYA 2008-02156 
and I. d. Gregorio by AYA 2008-06189 (Spanish grants). C. Pinte acknowledges 
funding from the EC 7''' FP as a Marie Curie Intra-European Fellow (PIEF- 
GA-2008-220891). J.-C. Augereau. G. Duchene, J. Lebreton, C. Martin-Zaidi. 
F. Menard and C. Pinte acknowledge PNPS, CNES and ANR (contract ANR-07- 
BLAN-0221) for financial support. D.R. Ardila, S.D. Brittain, B. Danchi, C.A. 
Grady, I. Pascucci, B. Riaz, A. Roberge, G. Sandell and CD. Howards acknowl- 
edge NASA/JPL for funding support. PACS has been developed by a consortium 
of institutes led by MPE (Germany) and including UVIE (Austria); KUL, CSL, 
IMEC (Belgium); CEA, DAMP (France); MPIA (Germany); IFSI, OAP/AOT, 
OAA/CAISMI, LENS, SISSA (Italy); lAC (Spain). This development has been 
suppoiled by the funding agencies BMVIT (Austria), ESA-PRODEX (Belgium), 
CEA/CNES (France), DLR (Germany), ASI (Italy), and CICT/MCT (Spain). 



References 

Brittain, S.D., Simon, Th. , Najita, J.R., and Rettig, T.W. 2007, ApJ, 659, 685 
Cutri, R.M., Skrutskie, M.F, van Dyk, S. et al. 2003, The IRSA 2MASS All-Sky 

Point Source Catalog, NASA/IPAC Infrared Science Archive 
Dent, W.R.R, Greaves, J.S., Coulson, I.M. 2005, MNRAS 359, 663 
Dent, W.R.R et al. 2010, PASP 

Dominik, C, Dullemond, C.P, Waters, L.B.FM. et al. 2003, A&A 398, 607 
Draine, B. T 2003, ApJ, 598, 1017 

Grady, C. A., Schneider, G., Hamaguchi, K., et al. 2007, ApJ, 665, 1391 
Juhasz, A., Bouwman, J. Henning, Th. et al. ApJ, submitted 
Kurucz, R.L. 1993, CD-ROM No. 13. Cambridge, MA, SAO 
Mathews, G.S. et al., this issue 

Meeus, G., Waters, L.B.FM., Bouwman, J. et al. 2001, A&A 365, 476 
Meeus, G., Bouwman, J., Dominik, C, et al. 2002, A&A 392, 1039 
Meeus, G., Juhasz, A., Helming, Th. et al. 2009, A&A, 497, 379 
Ott, S. 2010, in ASP Conference Series, Astronomical Data Analysis Software 

and Systems XIX, Y Mizumoto, K.-I. Morita, and M. Ohishi, eds., in press 
Panic, O., Hogerheijde, M.R., Wilner, D., and Qi, C. 2008, A&A 491, 219 
Pilbratt, G. et al. 2010, this issue 

Pinte, C, Menard, F, Duchene, G., & Bastien, P 2006, A&A, 459, 797 
Pinte, C, Harries, T J., Min, M., et al. 2009, A&A, 498, 967 
Poglitsch et al. 2010, this issue 

Raman, A., Lisanti, M., Wilner, et al. 2006, AJ, 131, 2290 
SiciUa-Aguilar, A., Hartmaim, L. W., Watson, D. et al., ApJ 659, 1637 



G. Meeus et al.: Gas in the protoplanetary disc of HD 169142: HerscheVs view 



Sloan G.C., Keller, L.D., Forrest, W.J. et al. 2005, ApJ632, 956 

Sylvester, R.J., Skinner, C.J., Barlow, M.J. et al. 1996, MNRAS, 279, 915 

Woitke, P., Pinte, C, Tilling, I., et al. 2010, MNRAS, in press 

Woitke, P, Kamp, I. Thi, W.-F. et al. 2009, A&A 501, 383 

Zactiarias, N., Monet, D.G., Levine, S.E. et al. 2004, BAAS Vol. 36, p.l418 

Zubko, V.G., Mennella V., Colangeli L. & Bussoletti E. 1996, MNRAS 282, 1321 



Dep. de Fisica Teorica, Fac. de Ciencias, UAM Campus 

Cantoblanco, 28049 Madrid, Spain 

e-mail: gwendolyn.meeus@uam.es 

School of Physics, University of Exeter, Stocker Road, Exeter EX4 

4QL, United Kingdom 

Universite Joseph Fourier - Grenoble 1/CNRS, Laboratoire 

d'Astrophysique de Grenoble (LAOG) UMR 5571, BP 53, 38041 

Grenoble Cedex 09, France 

UK Astronomy Technology Centre, Royal Observatory, Edinburgh, 

Blackford Hill, Edinburgh EH9 3HJ, UK 

SUPA, Institute for Astronomy, University of Edinburgh, Royal 

Observatory Edinburgh, UK Institute for Astronomy, University 

of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 

3HJ, UK 

School of Physics & Astronomy, University of St. Andrews, North 

Haugh, St. Andrews KY16 9SS, UK 

LAEX, Depto. Astrofisica, Centro de Astrobiologia (INTA-CSIC), 

PO. Box 78, E-28691 Villanueva de la Cafiada, Spain 

Institute for Astronomy (IfA), University of Hawaii, 2680 

Woodlawn Dr., Honolulu, HI 96822, USA 

Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, 

Leuven, Belgium 

SOFIA-USRA, NASA Ames Research Center, MS 211-3, Molfett 

Field, CA 94035, USA 

NASA Goddard Space Flight Center, Exoplanets & Stellar 

Astrophysics lab, Code 667, Greenbelt, MD 20771, USA 

Astronomy Department, University of California, Berkeley CA 

94720-3411 USA 

Eureka Scientic, 2452 Delmer, Suite 100, Oakland, CA 96002, USA 

ALMA, Avda Apoquindo 3846, Piso 19, Edificio Alsacia, Las 

Condes, Santiago, Chile 

European Southern Observatory, Alonso de Cordova 3107, Vitacura, 

Santiago, Chile 

Kapteyn Astronomical Institute, P.O. Box 800, 9700 AV Groningen, 

The Netherlands 

Spanish Virtual Observatory, Unidad de Archive de Datos, Depto. 

Astrofisica, Centro de Astrobiologia (INTA-CSIC), PO. Box 78, E- 

28691 Villanueva de la Caiiada, Spain 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St., 

Cambridge, MA, USA 

NASA Herschel Science Center, California Institute of Technology, 

Pasadena, USA 

Calar Alto Observatory, Centro Astronomico Hispano-Aleman, c/ 

Jesus Durban Remon, 2-2, 04004 Almeria, Spain 

Department of Physics & Astronomy, Clemson University, Clemson 

SC 29634-0978 USA 

NASA Exoplanet Science Institute/Caltech 770 South Wilson 

Avenue, Mail Code: 100-22, Pasadena, CA USA 91 125 

Max Planck Institut fiir Astronomic, Konigstuhl 17, 69117 

Heidelberg, Germany; 

Department of Physics and Astronomy, Johns Hopkins University, 

Baltimore, MD 21218, USA 

Research and Scientific Support Department of ESA, ESTEC/SRE- 

C, Postbus 299, 2200 AG Noordwijk, The Netherlands 

Astrophysikalisches Institut und Universitatssternwarte, Friedrich- 

Schiller-Universitat, SchillergaBchen 2-3, 07745 Jena, Germany 

Department of Radio and Space Science, Chalmers University of 

Technology, Onsala Space Observatory, 439 92 Onsala, Sweden 

ESA-ESAC Gaia SOC, PO. Box 78, 28691 Villanueva de la Caiiada, 

Madrid, Spain 

Spitzer Science Center, California Institute of Technology, 1200 E 

California Blvd, 91 125 Pasadena, USA. 



Department of Astronomy, Graduate School of Science, Kyoto 

University, Kyoto 606-8502, Japan 

CEA/IRFU/SAP AIM UMR 7158, 91191 Gif-sur-Yvette, France 

European Southern Observatory, Karl-Schwarzschild-Strasse, 2, 

85748 Garching bei Miinchen, Germany. 

The Rutherford Appleton Laboratory, Chilton, Didcot, OX 11 OQL, 

UK 

Department of Physics & Astronomy, The Open University, Milton 

Keynes MK7 6AA, UK and The Rutherford Appleton Laboratory, 

Chilton, Didcot, OX 11 OQL, UK 



